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How much are 2d Yukawa models similar to the Gross-Neveu models?* 
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We present numerical evidence that the 2d Yukawa models with strong quartic selfcoupling of the scalar field 
have the same phase structure and are asymptotically free in the Yukawa coupling like the Gross-Neveu models. 



1. Yukawa models in 2 dimensions 

The 2d Yukawa models (Y 2 ) with chiral Z(2) 
or U(l) symmetries are natural extensions of the 
usual or chiral 2d Gross-Neveu (GN) models, re- 
spectively. Starting from the auxiliary field repre- 
sentation of the 4-fermion coupling, one can add 
both the kinetic term and a self-interaction of this 
field 4> into the GN action. On the lattice the Z(2) 
symmetric Y 2 action is then (we introduce Np /4 
"naive" Dirac fermion fields ip a ) 
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(i) 



The 4> 4 scalar selfcoupling has been chosen, out of 
many possibilities in 2d, for the sake of analogy 
with the 4d theories. For n = A = the Z(2) 
GN model is obtained if the Yukawa coupling y is 
related to the usual GN coupling g by y = V2g. 

By choosing the above hopping parameter k 
formulation of the scalar field sector the kinetic 
term can be turned on or off gradually, elucidating 
the smoothness of the transition from the auxil- 
iary to a dynamical field. The relationship 
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between the Yukawa coupling y used in contin- 
uum and y is singular at k = 0, however. Another 
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virtue of the hopping parameter formulation is 
that spin models with the Z(2) or U(l) symme- 
try (the Ising or the XY models, respectively) are 
easily recovered at y = and k finite by choosing 
A = oo. The Y 2 models thus interpolate between 
the GN and spin models. 

2. Expected scaling properties 

Motivated by the recent discussion of a rela- 
tionship between the Nambu-Jona-Lasinio type 
four-fermion theories and the Standard Model |l]- 
U we address here the question in which regions 
of the parameter space the Y 2 models still possess 
the most cherished properties of the GN models, 
namely the asymptotic freedom at y — > 0, the 
fermion mass generation and, in the case of the 
Z(2) model, the symmetry breaking In the 

GN models these properties are derived by means 
of the 1 /Np expansion. 

For Y 2 models this expansion is applicable for 
A = 0(1/^^) 0,^] and here the same results as 
for the GN models are found. In particular, as 
long as k < k c (X = 0) = 1/4, the fermion mass 
amp is expected to scale as 



amp cx exp 



1 a 2 m 1 



2/3 Z yi 
For A — we have 



m = 777. n =11 — 



k c (0) / a 2 n ' 2k ' 



(3) 



(4) 



with too being the scalar field mass at A = y = 0. 
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For large A the 1 / Np expansion is a priori not 
applicable, not to speak of the perturbation the- 
ory. Nevertheless, M.A. Stephanov suggested |(| 
that in this case the mean field (MF) theory could 
be a good guide. An observation of long range fer- 
romagnetic couplings in the effective scalar field 
theory at y > by E. Seiler (?]] supports the ap- 
plicability of the MF theory. The resulting ex- 
pectations are ||: 

(i) The Y 2 models at A = oo or large pos- 
sess at arbitrarily small y only the broken sym- 
metry (in the Z(2) case) or the spin-wave (in the 
U(l) case) phase, whereas the symmetric or vor- 
tex phase present in the y = case in the pure 
scalar theory is absent at y > 0. 

(ii) The fermion mass and the magnetization 
y((j)} scale according to eq. (||). The coefficient 
m 2 is for k > the squared scalar mass at y = 
and Z its wave function renormalization constant. 
In general, Zj m 2 could be replaced by the scalar 
field propagator at zero momentum, i.e. the sus- 
ceptibility. This should hold for any k < k c (A) 
where n c (X) is the line of critical points in the 
pure scalar theory at y = 0. 

Thus according to the MF theory the asymp- 
totic freedom and the other mentioned properties 
of the GN models might occur also for large A in 
the Y2 models. 

3. Methods of data analysis 

We have tried three methods: 
(i) The asymptotic scaling law, which could 
be used to fit the data, is as (||) multiplied by 
y-Pi/0o, Here [3 and 0i are the /3-function coef- 
ficients which in the GN limit have the perturba- 
tive values &HT 
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in the U(l) case. This method does not take the 
finite volume effects into account. Nevertheless, 
one can roughly determine the values of the co- 
efficient Pq and compare with the values (pi) and 



(g). The 2-loop contribution given by /3i is less 
important than the finite size effects. 

(ii) A more appropriate way of analyzing the 
data is provided by a modified gap equation 
on finite lattices, 
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where the sum is performed over the fermion mo- 
menta with one periodic and one antiperiodic 
boundary condition on the L 2 lattices. Thanks 
to the IR divergence it leads for small amp in the 
infinite volume limit to the scaling law (^). 

The coefficient (3q is considered as a free pa- 
rameter in order to allow for its possible devia- 
tion from the perturbative value. The parameter 
s takes into account the fact that the gap equa- 
tion determines the value of the mass gap with 
insufficient precision. In order to fit the data it is 
necessary to treat the mass gap as a free param- 
eter. 

It would be interesting to determine the mass 
gap and compare it with the recently obtained 
exact results for the GN models @,|||. With 
our present accuracy it could be only estimated 
to be roughly consistent with these results. 

We note that one could take the finite size ef- 
fects into account also beyond the leading 1/Np 
order [ fT3||il| ] . As we do not expect this expansion 
to be applicable for large A, we do not attempt 
such refinements. 

Once the parameters in the gap equation have 
been determined by a fit, the quantity h(n, A), 
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has been extracted. 

(in) We have tried to fit the data also by an 
alternative to the essential singularity at y = 0, 
namely by a hypothetical power law behavior 



m F = a(y - y c ) 



4. Results at A = 



(9) 



To gain experience we have first performed sim- 
ulations at A = in the interval — 0.1 < k < 
k c (A = 0) = 1/4. We have determined in both 
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models the y-dependence of airiF and of y(<j>) in 
the Z(2) model at fixed values of ft on lattices 
of size 16 2 - 64 2 . The following observations are 
useful for the study of the models at large A: 

(i) The data analysis by means of the asymp- 
totic law (||) both without and with the 2-loop 
correction is possible if those points at small y, 
which show finite size effects, are excluded (fig. 1). 
However, the obtained values of (3q change with 
lattice size, when lower values of amp can be 
taken into account on larger lattices, so that one 
probably does not see the true asymptotic be- 
haviour. 

(ii) The gap equation ^ can describe the 
data obtained on different lattices consistently by 
means of one set of parameter values, including 
/3o (fig. 1). The onset of finite size effects at small 
arriF, as well as the general tendency of the data 
at large y are well reproduced. This analysis is 
superior to that by means of 

(Hi) The fermion masses in both models and 
y{4>) in the Z(2) model behave in a very analo- 
gous way. The magnetization in the U(l) case is 
present on finite lattices, but it shows a signifi- 
cant size dependence consistent with its vanish- 
ing in the infinite volume. Thus we observe the 
dynamical fermion mass generation taking place 
in spite of the absence of spontaneous symmetry 
breaking in 2d [||. 

(iv) The K-dependence is consistent with the 
expected one, eq. (|J). This is demonstrated in 
fig. 2, where /i(k,0) is shown. The values of (3q 
obtained at different k are for both models con- 
sistent with the values (||) and (^) within 10 %. 

(v) The power law formula (^|) can fit the 
data well for any given lattice size. The parame- 
ters b and y c depend strongly on the lattice size, 
however. This makes an algebraic singularity less 
probable than the essential one, but a reliable ex- 
clusion of the former on the basis of data alone 
seems very difficult. So one should be cautious 
at large A, when the analytic information on the 
type of singularity is less reliable than at A = 0. 

5. Results at large A 

Going to large A we have simulated both mod- 
els at A = 0.5 and the U(l) model at A = oo in the 



interval —0.3 < ft < k c (A). The most important 
results are: 

(i) The y-dependence of the fermion mass, 
including the onset of the finite size effects, is 
described by the gap equation (0) nearly as well 
as in the A = cases. This is demonstrated in 
fig. 3 for the U(l) model at A = oo and K = 0on 
16 2 - 64 2 lattices. 

(ii) The power law behaviour @ is not ex- 
cluded but, similarly to A = 0, disfavored by the 
strong dependence of y c and b on the lattice size. 
However, the finite size scaling analysis based on 
the Ansatz (||) has yet to be done. 

(Hi) These facts lead us to the tentative con- 
clusion that at large A, including A = oo, only 
the broken symmetry (Z(2) model) or spin wave 
(U(l) model) phase is present at arbitrary small 

y- 

(iv) Except at ft = 0, we do not yet have 
an independent determination of m 2 JZ for large 
A. Therefore we cannot yet determine the val- 
ues of flo and give the results in form of the val- 
ues of the coefficient 1i(k, A). As an example we 
show h(n,oo) for the U(l) model in fig. 4 (here 
k c (oo) w 0.56). We hope to determine (3q at large 
A for various ft in the near future. At present our 
estimates indicate that its values are roughly con- 
sistent with (3 in the GN cases, eqs. (|]) and (||). 

Thus we have found some evidence that the 
Y2 models with Z(2) and U(l) chiral symmetries 
with a strong 4 interaction behave for k < k c (A) 
as the GN models of the same symmetry, in par- 
ticular they are asymptotically free. The tran- 
sition from the GN to the spin model universal- 
ity classes takes place probably at small y and 
«; w ft c (A). 
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Figure 2. The n dependence of the coefficient 
h(n, 0) in the Z(2) model at X — is consistent 
with eq. (Q) (dashed line). 
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Figure 1. The 1/y 2 dependence of y ((f) in the 
Z(2) GN model (k = X = 0) with fits by eq. (||). 
The straight line is a fit by eq. (^). 



Figure 3. The l/y 2 dependence of amp in the 
U(l) Y2 model with Nf — 16 at A = 00 with fits 
by eq. §}. 
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Figure 4. The n dependence of the coefficient 
h(n, 00) in the U(l) model at A = 00. 



